The origin of cores and density profiles of gaseous baryonic 

structures 
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ABSTRACT 

We study the origin of cores and density profiles of gaseous baryonic structures in cosmology. 
By treating the baryons as a viscous gas, we find that both spheres and disks are possible 
solutions. We find analytically that the density profiles have inner and outer solutions, which 
in general are different. For disks we identify a central core, with density profile p e i = constant, 

r~ 6 . In the 



and the outer profile pa 



For spherical structures we find the profile p s 



presence of a dominating central black hole we find the inner profile p ~ r~ 3 / 2 . When the mass 
is dominated by a dark matter component then the baryonic density profile will depend on the 
dark matter profile, and we point out how one can use this connection to infer the DM profile 
directly by observing the baryonic density profile. 

Subject headings: galaxies: clusters: general — galaxies: structure — methods: analytical 



1. Introduction 

The universe is full of large gaseous baryonic 
structures such as galaxies and clusters. These 
structures are observed to take on a variety of 
shapes ranging from disks to spherical configu- 
rations. The radial profiles of these structures 
have been observed to be rather non-trivial, e.g. 
changing from one power-law profile near the cen- 
tre to another power-law profile at large radii. 
Even though these baryonic configurations have 
been observed and studied for many years, there 
is surprisingly little theoretical guidance to under- 
stand the origin of such complexity. Astronomers 
have for years been using phenomenological pro- 
files, ones with differing inner and outer shape, 
however, no simple theoretical explanation for this 
structure exists. 

In recent years numerical analyses have im- 
proved enormously, and we can now simulate much 
of the structures we observe. Never the less, it is 
important to have a simple understanding of the 
underlying physics, which can be obtained most 
easily through analytical studies of the basic equa- 
tions. 

In this letter we attempt an analytical treat- 



ment of the basic equations, asking which den- 
sity profiles are expected of purely gaseous bary- 
onic structures, and which profiles lead to stable 
baryonic structures within dark matter halos. We 
address these questions by considering a fluid ap- 
proach whereby we analyse asymptotic stable so- 
lutions of the Navier-Stokes equations. This ap- 
proach is new in astrophysics and surprisingly sim- 
ple. Our main results are, that the density profile 
of gaseous baryonic structures is fairly complex, 
namely that the radial density profile may have 
different slopes in the inner and outer regions. 
This is an issue which has been much discussed 
for Dark Matter profiles, but which has not been 
understood previously for baryonic structures. We 
find that spheres and disks are the only possible 
solutions. From an observational point of view 
it seems obvious that these are possible solutions, 
and while the creation of disks is understood theo- 
retically, the fact that no other solutions are stable 
is non-trivial. As a byproduct of our analysis we 
show how one can infer the DM density profile 
purely by observing the baryonic density profile. 
We also analyse how a central massive black hole 
influences the solutions. 



1 



2. Solving Navier-Stokes equations 

The behaviour of any collisional gas or fluid 
is fully determined by the Navier-Stokes (N-S) 
equations, which are 3 hydrodynamical equations 
for the velocity vector, and a continuity equa- 
tion. In the following we will make several as- 
sumptions, which certainly limits the applicabil- 
ity of our results, however, the solutions we find 
should be valid for the description of general prop- 
erties of structures such as gaseous baryonic galax- 
ies and haloes. The main assumption is that there 
are sufficient collisions while the structure forms. 
Thus if a gaseous baryonic structure has formed 
and subsequently all the baryons condense into 
stars, then the star density profile should still 
approximately follow the original density profile. 
Baryons often have sufficient collisions to equili- 
brate, e.g. in a typical intra-cluster gas the equi- 
libration timescale is about 10 7 yrs, with mean 
free path of tens of kpc compared to radii of few 
Mpc. We do not expect that our findings should 
apply to dark matter structures, since dark mat- 
ter presumably does not have sufficient collisions 
to ensure the validity of the N-S equations. 

In spherical coordinates the N-S equations de- 
scribe the velocity vector (v r , v@, v^). Here r is 
the radial coordinate, is the angle in the xy- 
plane, and cf> is the angle from the z-axis. The form 
of the equations is well known, see e.g. Landau 
& Lifshitz (1987), and contains time-derivatives, 
pressure terms, viscosity terms and a gravitational 
term. 

Our first 2 assumption are, that the gas 
has reached a stable configuration (no time- 
derivative), and that it has picked out an ori- 
entation in space, in such a way that all the gas 
is moving only in the 0-direction. Thus we have 
v r = (no contraction or expansion), and = 0. 
Here one must keep in mind, that by considering 
the N-S equations we are taking a fluid approach, 
which implies that we are following a fluid ele- 
ment, and this basically corresponds to averaging 
over all the particles moving through the fluid el- 
ement. For the 0- velocity we consider the general 
form 



we = v a 



(1) 



where a and \ are constants to be determined, v a 
and r a are unknown constants, with the physical 



interpretation that r a is a characteristic radius, 
and v a is the velocity of the fluid element at that 
radius. 

2.1. The O-equation 

The i>e-equation becomes very simple with the 
assumed form of the velocities 



= i/ 



V 2 w e - 



r 2 sm 2 c 



(2) 



where V 2 is the scalar Laplacian, and v is the 
kinematic viscosity. For now all that matters is 
the existence of a non-zero viscosity, so the ab- 
solute magnitude (and even radial dependence) is 
not important for the results . When we use the 
form for vq in cq. (1), then equation (2) has 4 
solutions 



X = +l 

x = -i 



a = 1,-2, 
a = 0, -1 . 



(3) 
(4) 



The solution with \ = +1 is exactly what one 
should expect for spherical symmetry, and we will 
refer to these solutions as the 'spherical solutions '. 
The solutions with \ = —\ indicate unstable solu- 
tions everywhere except in the flat cylinder where 
sin0 = 1, and we will refer to these as 'disk solu- 
tions'. We thus see that both spherical and disk 
solutions exist, and that they have different ro- 
tational structure (different a). Now, looking at 
eq. (1), it is clear, that the solutions with negative 
a are divergent for r — > 0, and we will therefore 
refer to those solutions as 'outer solutions ', and 
similarly, the solutions with positive a are incon- 
sistent for large radii, and we will refer to those as 
'inner solutions'. 

There are several obvious extensions of the 
method described here. One could in particular 
allow for a general O-dependence in tie, a non- 
zero radial velocity, or an exponential radial de- 
pendence of velocities and density profile 2 . 



1 Certainly baryons have a non-zero viscosity, however the 
viscosity could also arise from turbulence in which case it 
could be very large, f+urb ~ ^ v (Landau & Lifshitz 1987), 
where the dimension I is the size of the turbulent eddies, 
and At) is the velocity dispersion. 

2 Allowing a © dependence in vq as vq ~ (sinn0)' T , with n 
the number of spiral arms, gives <r = 0, 1 for all n. Allowing 
for a non-zero radial velocity, v r ~ r T , gives r = 1, —2. 
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Already from this one sees that a spherical so- 
lution, in eq. (3), must have different radial struc- 
ture for the inner and outer solutions (because the 
only solutions are a — 1, —2), and hence one may 
expect to find different density profiles in the cen- 
tral and outer regions. We point out that such 
phenomenon of simultaneous existence of two flow 
patterns is rather common in hydrodynamics, the 
simplest may be the hydraulic jump which is ob- 
served several centimetre large circular ring 
in any kitchen sink, when the flowing water goes 
from a 1/r profile to a constant, see Hansen et 
al. (1997); Watanabe et al. (2003). We emphasize 
that a non-zero viscosity appears to be a necessary 
condition for the existence of these specific solu- 
tions, even though a fundamental understanding 
of how the microscopic physics (viscosity) can de- 
termine the macroscopic properties (the general 
flow pattern) is still generally missing in hydrody- 
namics. 

2.2. The r-equation 

From the O-equation, we identified the general 
flow pattern, and we noticed the possibility that 
one may have different velocity-flows in the inner 
and outer region. We will now use the v r -equation 
to try to extract the asymptotic radial density pro- 
files. Also the r-equation is very simple 

v% ldP M(r)G 
r p or r z 

where p is the radially dependent density, P is 
the pressure, G is the gravitational constant, and 
M(r) is the mass within the radius r. We assume 
that the pressure and density are related through 
P = P a (p/ pa) 1 , where P a and p a are the un- 
known pressure and density at r a . We assume the 
gas is monatomic with 7 = 5/3. Let us consider 
densities of the form 

p(r)=p a (j-y, (6) 

such that the parameter (3 determines the density 
profile, ft is worth emphasising that it is exactly 
this (3 which we are trying to find. 

Let us study the radial dependence of the 3 
terms in eq. (5). Using v B in eq. (1) the first (ki- 
netic) term of eq. (5) goes like Vg/r ~ r 2a_1 . The 
pressure gradient term goes like 1/p ■ dP/dr ~ 



r 5 ' 3-1 , where we have used 5 = 7 — f = 2/3. 
The last (gravitational) term including M(r), de- 
pends on the given system we are considering. If 
the mass is dominated by a point gravitational 
source (e.g. a central black hole), then it goes like 
M(r)G/r 2 ~ r~ 2 . If the mass is dominated by the 
matter density, then it goes like M(r) ~ J p(r)dV , 
with dV the volume element. For spherical solu- 
tions this gravitational term thus goes like r^ +1 
with (3 from eq. (6), and for disk solutions it 
goes like r /3 '. This gravitational term has the cor- 
rect form for spherical distributions (and point 
sources) , but is only an approximation for the pure 
disk case. Technically speaking the mass is log- 
arithmic divergent for spherical structures with 
(3 = —3, however, the formula M ~ r@ +3 holds 
for any f3 arbitrarily close to —3, and furthermore 
in a real situation there would be an outer cut-off. 
To be explicit, we are looking for solutions to an 
equation of the form 

r 2a-l =J .Sf}-l +r . Kj ( 7) 

where k = (3 + 1, (3, —2 for spherical, disk and BH 
matter dominance respectively. When we use the 
word 'solve' in the following, we are really just 
using the standard method of divergence cancel- 
lation, in the sense that the most divergent terms 
must cancel with each other. The optimal case 
is naturally that all divergences disappear, a case 
which we will refer to as 'good'. From our simple 
analysis the transition radius, r Q , which separates 
the inner from outer region, is not uniquely deter- 
mined. We only find approximate disk relations 
like r a ~ v^/Gp a . A full study including the co- 
efficient is significantly more involved, and we will 
leave that for a later analysis 3 . 

We want to solve eq. (5) for the profile param- 
eter (3, however, there are 4 situations to consider 
(the 4 different a from eqs. (3,4)), and for each 
case we can choose which is the dominating mass 
contribution, either from a sphere, from a disk or 
from a black hole (BH). The discussion in sections 
2.3 and 2.4 is somewhat technical, and the reader 
is encouraged to look at figs. 1 and 2, which sum- 
marize the main findings of these sections. 

3 We note that when the Vq term is negligible this eq. (5) is 
exactly the hydrostatic equilibrium, which is often written 
as M(r) ~ rT^y/3. Furthermore, in the same notation one 
finds that the velocity dispersion, cr^, has a radial depen- 
dence, dlncr^ I dlnr = 8/3. 
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2.3. Disk 

Let us first consider the a = case, which can 
be both an inner and outer solution. For the inner 
solution we find two density profiles, (3 = 0,-3. 
The (3 = solution is more natural, because the 
first two divergent terms in eq. (5) (or similarly 
cq. (7)) can cancel, and the last term is non- 
divergent for r — > 0, whereas for the (3 = — 3 solu- 
tion, the two most divergent terms cancel, but the 
remaining term is still slightly divergent towards 
the centre. With this kind of argument one can 
divide all the solutions into 'good', and 'reason- 
able' solutions, and we will emphasise which are 
good solutions. If the mass is dominated by a BH, 
then the solution is (3 = —3/2. If the disk is dom- 
inated by a spherical distribution (cither from a 
large fraction of the baryons, or from an unknown 
dark matter component) with profile (3 S , then the 
solution is (3 = 3/2 (f3 s + 2). In conclusion we 
see, that the only good solution for the inner disk 
is (3 = 0. In principle the a = case can also 
be an outer solution. In that case the solution is 
(3 = — 1. If the mass is dominated by a BH, then 
[3 = (good), and if dominated by a spherical 
distribution then (3 < 0. 



Disk solutions 
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Fig. 1. — A summary of disk solutions. The full 
lines (self, (3 = 0, —3) are where the mass is dom- 
inated by the disk itself. The dotted line (BH, 
(3 = 0) is for black hole domination, and the 
dashed line (sphere, (3 = —3) is for matter domi- 
nation by a sphere with profile , (3 S = —3. 

The a = — 1 (which is an outer solution) gives 
(3 = — 3 (good). If a BH dominates then (3 = 
—3/2, and if another spherical component domi- 



nates (with profile (3 S ) then (3 = 3/2 ((3 s +2), where 
the case (3 = — 3 is the only good solution. We 
summarize the good disk solutions in Figure 1. 

To conclude the disk solutions, when the mass is 
dominated by the matter itself, then there is only 1 
good configuration (we emphasise that by good we 
mean optimal removal of divergences), which has 
the inner slope of (3 = and outer slope of (3 = 
—3. Following the similarity with the hydraulic 
jump in the kitchen sink, we will refer to such 
configuration as a galactic hydraulic drop. 

2.4. Sphere 

For the sphere it turns out that almost all the 
solutions are good, in the sense that all the di- 
vergences cancel in a simple manner. The a = 1 
is the inner solution, for which we find (3 = 3, —6. 
The (3 = 3 is obviously strange, and probably non- 
physical (a positive (3 would lead to a wrong sign 
in front of the pressure gradient term). With BH 
dominance one has (3 = —3/2, and if dominated 
by a disk (which could either be from a fraction of 
the same gas, or from another particle type) with 
profile (3d one finds (3 = 3/2 ((3d + 1). If dominated 
by another spherical distribution (which probably 
should arise from another particle type) with pro- 
file (3 S , then we find (3 = 3/2 ((3 S + 2). 

The outer solution, a = — 2 gives (3 = —6. If a 
BH dominates then (3 = —3/2, and again if a disk 
(or sphere) dominates then (3 = 3/2 ((3d + 1) (or 
3/2 ((3 S + 2)). We summarize the good spherical 
solutions in Figure 2. 

We note that if the mass is disk dominated with 
disk-profile (3d = —3, then the sphere happens to 
also get the profile (3 S = —3. Thus one can envis- 
age a significant fraction of the matter collapsing 
to a disk, which very well could take the profile 
(3d = —3 (as shown in section 2.3), and this would 
force the remaining matter, which is in a spherical 
configuration, to take the density profile (3 S = —3. 
We also point out that if the mass is dominated 
by a dark matter sphere, then the baryon pro- 
file becomes (3 = —3/2 for a dark matter profile 
(3 S = —3 as suggested by dark matter N-body sim- 
ulations. If instead the dominating DM profile is 
(3 S > —2, this would seemingly imply a positive 
baryonic slope, which is difficult to interpret. We 
feel that this limited applicability warns, that a 
more general calculation may provide a different 
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Spherical solutions 
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Fig. 2. — A summary of spherical solutions. The 
full lines (self, (3 = —6) arc where the mass is dom- 
inated by the sphere itself, and the dotted lines 
(BH, (3 = —3/2) are for black hole domination. 
The dashed lines (sphere, (3 — —3/2) are for mat- 
ter domination by another sphere with profile -3, 
and the dot-dashed lines (disk, (3 = — 3) are for 
domination of a disk with profile -3. 

connection between the DM and gaseous baryonic 
profiles 

3. Discussion 

An interesting possibility now appears, namely 
that one can use our results to infer the dark 
matter profiles from observations of the baryonic 
profile. This method can be used quite gener- 
ally to infer the DM distribution, and is there- 
fore complementary to other methods such as lens- 
ing observations. Let us say we have observed a 
baryonic sphere with profile /Sbaryon, and we know 
that the mass is dark matter dominated. Under 
the assumption that the DM is spherical we have 
Aaryon = 3/2 ((3 DM + 2). Thus, if we observe e.g. 
/Sbaryon = —1-5, then we know, that the DM has 
profile (3dm = —3. Determining the baryonic den- 
sity profile directly, i.e. independent of dynamics, 
can be accomplished in several ways: 

1) X-rays: observations of the luminosity, L x , 
in various bands and in different radial bins gives 
the radial electron density of the plasma. Here the 
main concern is the validity of hydrostatic equi- 
librium and disentangling any cooling flow in the 
centre of the cluster. Relaxed clusters with no evi- 
dence of cooling flows do exist, for example A2029 



(see Figure 1 of (Lewis et al. 2003)). The gas in 
this cluster shows an inner profile Sbaryon = —0.55 
and an outer /3baryon = —1-62 (note our differing 
definition of [3 here). The outer value is certainly 
DM dominated and our analysis implies that (3dm 
is very close to —3, which is expected from CDM 
simulations (Navarro ct al. 1996; Moore et al. 
1999). We do not attempt to deduce the inner 
DM slope here, because one cannot be certain that 
the mass is DM dominated at such small radii 
where the baryonic slope should reach its asymp- 
totic value. Again we refer to the warning in the 
end of section 2.4. 

2) The Sunyaev-Zeldovich effect is in principle 
a direct measure of the plasma column density. 
The angular resolution required for our analysis 
currently limits the use of this technique, however, 
future multi-frequency observations will determine 
independently both the temperature and number 
density profile of distant clusters purely through 
the S-Z effect (Hansen et al. 2002; Aghanim et al. 
2003). The SZ-effect can thus be used to measure 
the electron density profile to large radii (SZ effect 
is proportional to n e , whereas X-ray is propor- 
tional to n 2 e ) and at large redshift (SZ is redshift 
independent). 

3) Surface brightness: both from stellar light 
and radio observations of HI and molecular gas, 
one can in principle determine the baryon profile. 
An example is for M33 by (Corbelli 2003). 

The distribution near black holes has been con- 
sidered earlier. First by Peebles (1972) where en- 
ergy consideration lead to the distribution p ~ 
r~ 9 / 4 , which was refined in (Bahcall & Wolf 1976) 
who found p ~ r~ 7 / 4 . Later numerical simulations 
have shown (Young 1980) that the profile near the 
black hole should be p ~ r~ 3 / 2 , which is just what 
we find. For accreting black holes this profile may 
be different (Freitag & Benz 2002). 

4. Conclusion 

We make a first attempt to derive analytically 
the asymptotic density profiles of baryonic struc- 
ture in cosmology, which include galaxies, gaseous 
haloes and intra-clusters gas. We find that both 
disks and spherical solutions exist, and that gen- 
erally the inner and outer density profile may be 
different. Thus we supply theoretical support for 
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the use of phenomenological profiles like Watanabe, S., Putkaradze, V. & Bohr, T. 2003, J. 

Fluid Mech., 480, 233 

= r 1 (i +r y 2 ■ ( 8 ) Young, P. 1980, ApJ, 242, f 232 

For the disks we identify central cores with (3d = 
from p ~ r' 3 . For the outer region we find (3d — — 3. 
For spherical structures we identify both inner and 
outer profiles, which include (3 S = —6,-3,-3/2. 
Our resulting profiles only apply to gaseous bary- 
onic structures, but we point out a simple method 
whereby observations of the baryonic structure in 
principle allows one to deduce the dark matter 
density profile. 
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